High Temperature Universal Properties of Atomic Gases Near Feshbach Resonance 

with Non-Zero Orbital Angular Momentum 
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We show that the high temperature behavior of atomic gases near Feshbach Resonance with 
non-zero orbital angular momentum {£ > 0) belong to a universality class different from that of 
s-wave resonances. The universal interaction energy is 2(2^-1- 1) times larger than that of the s-wave 
when approaching the resonance from the atomic side, but is essentially zero on the molecular side; 
contrary to s-wave resonances where interaction energies on both sides are the same except for a sign 
change. The measurement of these universality properties should be feasible in current experiments. 



In the last eighteen months, Feshbach resonance 
has been used with great success to achieve molecular 
condensates P and fermion superfluids0]. At the same 
time, the physics of Feshbach resonance poses a chal- 
lenging many-body problem due to its non-perturbative 
nature. Through Feshbach resonance, a pair of atoms 
can be converted into molecules. Just before and af- 
ter a bound pair is formed, the scattering length as di- 
verges, changing from positive infinity to negative infin- 
ity across the resonance. These divergences render the 
usual perturbative scheme in terms of the gas parameter 
n^/^Os inapplicable (where n is the density), and is the 
source of difficulty in theoretical treatment. At the same 
time, it implies that the system can exhibit "universal be- 
havior" , provided there are no other anomalously large 
length scales in the system. The reason is that the di- 
verging scattering length must disappear from the phys- 
ical properties of the system. If all other length scales 
are smaller than the inter-particle spacing and the 

thermal wavelength A, the thermodynamic functions of 
the system can only depend on these lengths and not on 
any microscopic properties, and is in this sense univer- 
sal. Such universal behavior has indeed been observed in 
many recent experiments 

At present, most experiments are on s-wave Feshbach 
resonance. However, there are many other resonances 
with non-zero orbital angular momenta {£ > 0). Very 
recently, Salomon's group at ENS 5] has reported a re- 
versible production of molecules in Fermi gas of ^Li across 
a p-wave resonance, paving the way for condensation of 
p-wave molecules and realization of p-wave fermion su- 
perfluid in the future. There is also a very exhaustive 
recent study of the Feshbach resonances of the Bose gas 
^^•^CsQ. Resonances upto 1 = 4 have been observed. 
Many of the molecules appear to have long lifetimes!^. 
It is natural to ask how quantum gases near £ > reso- 
nances differ from those of s-wave resonances. Will there 
be universal behavior near resonance? And if there is, do 
they belong to the same class as s-resonances. What is 
the nature of the £ > fermion superfluids in the strongly 
interacting region, and what is the nature of the molecu- 
lar condensates? In this paper, we shall address the issue 
of universality near £ > resonance. The studies of the 
£ > fermion superfluids and molecular condensates are 
quite involved and will be discussed separately. 



Since the physics at resonance is non-perturbative, it is 
useful to have exact results that can be used as a guide for 
theoretical approximations. Such exact results are pos- 
sible in the high temperature regime, where the grand 
potential can be expanded in powers of the fugacityQ. 
Note that the conceptual problem related to the diverg- 
ing scattering length does not disappear at high temper- 
atures. For example, it is well known that the interaction 
energy density for s-wave scattering far from resonance is 
gn^ at all temperatures, where g — AttH^ as / M , and M is 
the mass of the atom. This expression can not persist at 
resonance because of diverging g. In the case of s-wave, 
one of us has shown recently that the interaction energy 
density at high temperatures will change from gn? to the 
universal value (3nfc_BT/2)(nA'^/2^/^) as one approaches 
the resonance 0- Here, we shall perform similar exact 
calculations for £ > resonance. We have in mind a 
two component Fermi gas with identical number in each 
component, n-f = = n/2. Our exact results should be 
useful for future experiments. 

Our calculation reveals several remarkable features: 
(A) At high temperatures or low densities, the inter- 
action energy densities Cj^^. for all £ > resonances have 
identical behavior across resonance, apart from a trivial 
degeneracy factor {2£ + 1). They are, however, different 

from that of s-wave. (B) If we denote e|^^ and e|^^ the 
interaction energy density when approaching the reso- 
nance from the atomic or molecular side, (defined as the 
side of resonance where bound state is absent or present 
respectively), then for ^ = 0, e|°j and e|^'' are antisym- 
metric about the resonance. Such symmetry is lost for all 
£ > resonances. There, ej™"* is essentialy zero, whereas 

e[°j is large and negative close to resonance, though it 
decreases rapidly away from resonance. (C) For s-wave 
scattering, the second virial coefficient 62 and hence e\^^ 
reach universal values 1/2 and — (3nfcBT/2)(nA^/2^/^) 

respectively at resonance. For ^ > 0, 62 and e|^j are 
twice as big, in addition to a degeneracy factor (2£ -I- 1). 
These differences are due to the presence of the centrifu- 
gal barrier in the £ > scattering channels. These results 
are derived below. 

(I) Second virial coefRcient and interaction en- 
ergy density: Let us recall the well known high tem- 
perature expansion of the grand partition function Z = 
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FIG. 1: s-wave phase shifts for a square weU potential with 
width 6 and various depths. Different curves starting from 
bottom to top correspond to Oo/b = —5, —20, —100, ±oo 
(dotted-dash line), 100, 20, 5. The effective range at reso- 
nance is ra = h. A typical thermal wave length within the 
universal range of temperatures is also indicated. 
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FIG. 2: Similar to Fig. 1, only for the p-wave phase shifts. 
The successive curves starting from bottom to top are ai/& 
= -5, -20, -100, ±oo (dotted-dash line), 20, 5, 1. For the 
square well potential the p-wave effective range at resonance 
is ri = —36. 



'Pj.g-(ff-/iiV)/feBT powers of the fugacity z = e^l^^T ^ 
It is shown by Beth and Ulenbeck in 19370 that Z ~ 
Z(°)+2^/2 62 + 0(2^) + ... where "o" denotes quan- 

tities for the non-interacting system, V is the volume 
of the system, and 62 is the second virial coefficient de- 
fined as 62 = (e-'^^^" - e-Z^I^^"!'"') , where S^'^ and 

\E^^^''^ are the energy eigenvalues of a two-particle sys- 
tem with and without interaction. Separating the bound 
states from scattering states, Beth and Ulenbeck showed 
that 62 can be expressed asjSj 
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where E\, is the energy of the two-body bound state, 
biik) is the phase shift of the ^-th partial wave, and 

A = \l 2-kTi' I {MIibT) is the thermal wavelength. From 
the fugacity expansion, it is straightforward to show (7(| 
that the energy density can be expanded in the small pa- 
rameter riA'^, and e{T,n) = e]^in{n,T) + eint(n, T), where 
ekin(n, T) = 3n|BT ± + ...^ ^ is the kinetic energy 
of an ideal gas, and 



Cint = [nX ) 



V2 9^2 
V2 3 dT 



(2) 



is the interaction energy density. The nX^ term in ekin is 
the well known statistical correction, with -|- and — sign 
for bosons and fermions respectively. 

(II) Effect of centrifugal barrier on the bind- 
ing energy and the phase shifts: To understand 
the general features of the interaction energy near res- 
onance, let us recall the relation between scattering am- 
plitude f{k,d) and the phase shift Se{k), where 6 is the 
angle between incident and scattered waves. Resolving 
into partial waves, /(fc, 9) = ^g{2i+l)Pf {9)f£{k), where 



Pi{0) are the Legendre polynomials, it is well known that 
fi{k) is related to the phase shift S{k) of the £-th partial 
wave as fi{k) — k~^{cotS£ — For kb « 1, where 

b is the range of the scattering potential, 5e{k) has an 
expansion 

(kb)-^'^ f 1 , 



cot 5(,{k) 



k 



kb « 1, 



(3) 

where ai and rg are the scattering length and effective 
range of the i-th partial wave respectively. The low en- 
ergy scattering amplitude can then be written as|lCl| 

(kb)^'- 



h{k) = 



-l/ae + rek'^/2-ik{kb) 



21 ' 



fc6 « 1, (4) 



It is also well known that by analytically continuing k 
in the complex plane, if /^(fc) has a singularity on the 
positive imaginary axis, k = in, k > 0, then the system 
has a bound state with energy 

The factor (kb)'^^ in eq.Q is due to the presence 
of a centrifugal barrier. It has dramatic effects on 
the bound state energy and phase shift that appear in 
eq.(^. For s-wave scattering, as fc ^ 0, eq.|0J becomes 
/o(fc) ^ l/(^'^o ^ ~ *^)- The dominant k dependence is 
the imaginary part in the denominator. When k is an- 
alytically continued to the pure imaginary axis, k — in, 
with K > 0, /o(fc) has a singularity (and hence a bound 
state) only when oq > 0; in which case k = I/oq and 
E\j — ~h^/{Maf^). In contrast, for ^ > 0, the low en- 
ergy behavior of fi in eq.(0J) is dominated by the real 
part of the denominator. A bound state occurs when 
fc^ = —K^ with K — ■\/2/(— a^r^), which is possible only 
when a^Tf, < 0. The binding energy is 
Unlike the s-wave case, Ei, is linear in l/ag. In addi- 
tion, it depends on effective range rg which is typically 
of atomic scale. 

Turning to the phase shifts, eq.(|3J) implies that 



d5t{k) _ [kbf^ [-{2i+l)a~^ + (2£- l)r^fcV2] 



(5) 
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FIG. 3: 62 for £^0. For ^Li and b = lOnm ^ 200aB the cor- 
responding temperatures are dashed, sohd and dotted hnes 
50/iK, 12^K and 2fj,K. The deviation of 62 from ±1/2 at res- 
onance is due to effective range corrections. 



FIG. 4: eint for £ = 0. The same temperatures as in Fig. 3. 
fcF is related to density as k% = 3n^n, which we choose to be 



When £ ~ 0, the denominator is k"^ + (og ^ — roA;^/2)^. 
Near resonance |ao| >> tq and one can ignore ro to sec- 
ond order in k, and hence 



dSojk) 
dk 



-ao 



1 + alk^ 



+ O (rofc2) 



(6) 



which is a half Lorenzian in the range of positive k. Note 
that this result is valid on both sides (molecular and 
atomic) of the resonance. In contrast, for £ > 0, the 
term (kb)'^^ is negligible in the denominator of eq.©, so 
the effective range can not be ignored. This makes 
dSe/dk very different on the two sides of the resonance. 
On the atomic side where bound states are absent, we 
have a^ri > 0. Near k = kc = ■\/2/ (airi), we have 



d5i{k)_ 
dk 



marked difference in behavior lead to the differences in 
dSg/dk shown in eq.© and {T)). 

(Ill) Final Results: Let us first explain the reason 
for the emergence of universal behavior near resonance 
before deriving the accurate formulas. The first term in 
eq. is certainly universal near resonance since i?b — s- 
as the bound state emerges. The second term in eq.([3) 
is the contribution due to scattering states. This integral 
is cut off by 1/A, and is roughly given by 

,«..,a+i)r"iiMfi.fi«)[,,,,-.)_*(o)] 

Jo '^'^ 

.(8) 

Near resonance, the s-wave phase shift (5o(^) quickly rises 
from zero to 7r/2 and stays basically flat within the region 
« A: « \ where < cot(5o « 1. (See eq.® 



-1/2 



(fc - kcY + Tiikcf 



air I > 0, (7) and also Fig. 1). We then have 



where r<>(fc) = \{kb)^'^ /rp]. Eq.0 is a very narrow Lozen- 
zian because Tt{kJ)/kc = \{2h/ atY-^''^{b/ri.Y+^/^\ « 1. 
Thus, unlike the s-wave case where dSo{k)/dk is a half 
Lorenzian in the positive fc-axis, di5^(fc)/dfc is a /mH Loren- 
zian in the same range of k for £ > 0. On the molecular 
side where a^r^ < 0, dS{k)/dk has no resonance struc- 
ture. Since it is proportional to cx {kbY^, Si is exceed- 
ingly small. 

The differences between 60 and 5e are shown explicitly 
in Fig. 1 and Fig. 2 (for £ — 1). Let us consider both 
cases very close to resonance on the atomic side. For 
£ — 0, Sq rises quickly from to 7r/2 as k increases. It 
stays close to 7r/2 within the region Oq ^ « k « r^^ , 
where cot i5o satisfies < cot Sq « 1 (see eq.©). It then 
begins to decrease as k increases. (Note that all Si{k) 
must vanish as k — > 00, since the effect of scattering 
potential becomes unimportant at very high energies). 
For ^ > 0, however, 6e{k) remains close to zero until k 
approaches kc = ■\/2/(a^r^). It rises quickly from to tt 
as k passes t hrough k c- It stays close to tt — 0+ within 
the region yJYfia^r^ « k « b'^\re/b\^/'-^^~'^^ where 
cotSe « —1, and then decreases as k increases. These 



1/2 for To « A « n 



-1/3 



(9) 



The condition Xn^^'^ « 1 is automatically enforced be- 
cause of the low fugacity,7j. Eq.® shows that univer- 
sal behavior only emerges when the temperature is suffi- 
ciently high to achieve low fugacity limit but sufficiently 
low so that the thermal wavelength is still much larger 
than the effective range, which is typically the range of 
the potential. In fact, at very high temperatures where 
A ^ 0, 62 vanishes for all £ since S{X~^) S{oo) = 0. 
However, such high temperature range is not of interests 
in current experiments. 

For £ > 0, 5i{k) rises rapidly from to tt as fc passes 
through kc = y^2/{aere). It stays close to tt within 
the region « k « b^'^{re/by/'-^^^'^\ where 

cot Si « —1. (See also fiure 2). This means that 

62 » (2^+1) for 6(6/r£)i/(2^-i) « A « (10) 

Having explained the emergence of universality near 
resonance, we now derive the precise formula for interac- 
tion energy for all region of scattering length. The case 



4 





b«/3 

1 

\ ^\ 


A,=10b 

— A,=20b 

— A,=50b 


0.5- 


\ ^\ 

\ 




0.01 


0.01 


0.02 0.03 



-b/a 



FIG. 5: 62 for 1=1. The dashed, solid, and dotted lines 
corresponds to T = 50, 12, 2/LtK. The deviations of 62 from 1 
on the atomic side and on the molecular are due to effective 
range corrections. 



where erf (cc) is the error function and x = A/\/27rao. The 
behaviors of 62 and eint calculated from eq. ^ and Q are 
shown in figures 3 and 4. The results in these figures are 
exact calculations for the square well potential and are 
found to be indistinguishable from eq. Hill) when effective 
range correction is included. The antisymmetry of 62 
and Cint on different sides of the resonance is due to the 
antisymmetry of eq.(|SJ in I/oq. From eq.JQ) and jSJ, one 
can easily show that at resonance, 62 = ±1/2 — 
which confirms eq.®. The corresponding universal value 
for the energy density is eint = {3nkBT/2){nX^ /2^/^) + 
O{ro/X). 

For £ > 0, eq.© and imply that on the atomic side 



62 = {2l+l)e~ 



""ff erfc 
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FIG. 6: tint for £ = 1. Same parameters as in Fig. 4. Note 
that the molecular branch (top left curves) is close to 
and the atomic branch at resonance is larger than that of the 
s-wave. 



of s-wave resonance has already been discussed in ref . ■ 
Substituting eq.® into eq.Q), we have for £ = 



&2 = 5]e" 



% sgn(ao) 



1 — erf (a;) 



O 



(11) 



where erfc(x) = 1 — erf(x), and (...) means correction 
of the order of [j]^*^ The behaviors of 62 

and eint for £ = 1 are shown in figures 5 and 6. They 
are exact calculations for a square well and are in ex- 
cellent agreement with ea. (ll2|l . Approaching the res- 
onance from the atomic side, we see from eq.lO that 
&2 = {2£+l) (1 - {2£ - l)!!7r'^23/2 A [ &] ^^-^^ ^ which ver- 
ifies ea. ((TU|) . The interaction energy density is eint = 
2(2f + l)(3nfcsT/2)(nAV23/2) +0 (A [^j^^" 

On the molecular side, eint is zero to the zeroth or- 
der in ro/ A because of the exceedingly small phase shifts 
(~ k'^^'^^). Far from resonance, 5i{k) — —aik{kb)'^^, 
we have for all £ and on both sides of the resonance, 
42 = gen^n^ (^f ' ^^(2£-f 1)!! ([1 - £] - ^^£i2£ + 3)), 
where gi — AttTi a^/M^ and = ni = n/2. This shows 
that away from resonance, the interaction energy in the 
£ > Q channels can not be described by a single length 
scale ai as in the case of s-wave scattering. We have thus 
established the results (A) to (C) in the Introduction. 
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